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Spontaneous  symmetry  breaking  can  lead  to  the  formation  of  time  crystals,  as  well  as  spatial  crystals. 

Here  we  propose  a  space-time  crystal  of  trapped  ions  and  a  method  to  realize  it  experimentally  by 
confining  ions  in  a  ring-shaped  trapping  potential  with  a  static  magnetic  field.  The  ions  spontaneously 
form  a  spatial  ring  crystal  due  to  Coulomb  repulsion.  This  ion  crystal  can  rotate  persistently  at  the  lowest 
quantum  energy  state  in  magnetic  fields  with  fractional  fluxes.  The  persistent  rotation  of  trapped  ions 
produces  the  temporal  order,  leading  to  the  formation  of  a  space-time  crystal.  We  show  that  these  space- 
time  crystals  are  robust  for  direct  experimental  observation.  We  also  study  the  effects  of  finite 
temperatures  on  the  persistent  rotation.  The  proposed  space-time  crystals  of  trapped  ions  provide  a 
new  dimension  for  exploring  many-body  physics  and  emerging  properties  of  matter. 

DOl:  10.1 103/PhysRevLett.l09.163001  PACS  numbers:  37.10.Ty,  03.65. Vf,  64.60.Bd,  64.70.Nd 


Symmetry  breaking  plays  profound  roles  in  many-body 
physics  and  particle  physics  [1,2].  The  spontaneous  break¬ 
ing  of  continuous  spatial  translation  symmetry  to  discrete 
spatial  translation  symmetry  leads  to  the  formation  of 
various  crystals  in  our  everyday  life.  Similarly,  the  sponta¬ 
neously  breaking  of  time  translational  symmetry  can  lead 
to  the  formation  of  a  time  crystal  [3,4].  Because  of  sym¬ 
metry  breaking,  the  effective  ground  state  of  a  system  can 
be  an  inhomogeneous  many-body  state  whose  probability 
amplitude  changes  periodically  in  time  [3],  having  a  tem¬ 
poral  symmetry  different  from  that  of  the  true  quantum 
ground  state  of  the  Hamiltonian.  This  is  similar  to  the 
symmetry  breaking  during  the  formation  of  a  spatial  crys¬ 
tal,  where  the  crystal  has  a  symmetry  different  from  that  of 
the  underlying  system  [5]. 

Intuitively,  if  a  spatially  ordered  system  rotates  persis¬ 
tently  in  the  lowest  energy  state,  the  system  will  reproduce 
itself  periodically  in  time,  forming  a  time  crystal  in  analog 
of  an  ordinary  crystal.  Such  a  system  looks  like  a  perpetual 
motion  machine  and  may  seem  implausible  in  the  first 
glance.  On  the  other  hand,  it  has  been  known  that  a 
superconductor  [6,7]  or  even  a  normal  metal  ring  [8-10] 
can  support  persistent  currents  in  its  quantum  ground  state 
under  proper  conditions.  However,  the  rotating  Cooper 
pairs  or  electrons  in  a  metal  are  not  quantum  time  crystals 
since  their  wave  functions  are  homogeneous  and  no  time 
translational  symmetry  is  broken.  A  soliton  model  that 
assumes  strong  attractive  interaction  between  particles 
with  the  same  type  of  charge  has  been  proposed  [3],  which 


is  difficult  to  be  realized  experimentally.  While  it  has  been 
proved  mathematically  that  time  crystals  can  exist  in  prin¬ 
ciple  [3,4],  it  was  not  clear  how  to  realize  and  observe  time 
crystals  experimentally. 

In  this  Letter,  we  propose  a  method  to  create  a  space- 
time  crystal  [Fig.  1(a)],  which  is  also  a  time  crystal,  with 
cold  ions  in  a  cylindrically  symmetric  trapping  potential. 
Different  from  electrons  in  conventional  materials,  ions 
trapped  in  vacuum  have  strong  Coulomb  repulsion  be¬ 
tween  each  other  and  have  internal  atomic  states.  The 
strong  Coulomb  repulsion  between  ions  enables  the  spon¬ 
taneous  breaking  of  the  spatial  translation  symmetry,  re¬ 
sulting  in  the  formation  of  a  spatial  order  that  can  be 
mapped  into  temporal  order.  The  internal  atomic  states  of 
ions  can  be  utilized  to  cool  the  ions  to  the  ground  state  as 
well  as  observing  their  persistent  rotation  directly  by  state- 
dependent  fluorescence.  Different  from  Wilczek’s  soliton 
model  [3],  our  trapped-ion  model  is  a  crystal  (in  space) 
even  when  the  magnetic  flux  is  zero. 

Trapped-ion  Coulomb  crystals  have  provided  unique 
opportunities  for  studying  many-body  phase  transitions 
[11-13]  and  quantum  information  science  [14,15].  The 
two  most  common  types  of  ion  traps  are  the  Penning  trap 
[16]  and  the  Paul  trap  [15].  The  magnetron  motion  of  an 
ion  in  the  Penning  trap  is  an  orbit  around  the  top  of  a 
potential  hill  that  is  not  suitable  for  studying  the  behavior 
of  the  ion  in  its  rotational  ground  state.  Here  we  propose  to 
use  a  combination  of  a  ring-shaped  trapping  potential  from 
a  variation  of  the  Paul  trap,  and  a  weak  static  magnetic 
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FIG.  1  (color  online).  Schematic  of  creating  a  space-time 
crystal,  (a)  A  possible  structure  of  a  space-time  crystal.  It  has 
periodic  structures  in  both  space  and  time.  The  particles  rotate  in 
one  direction  even  at  the  lowest  energy  state,  (b),  Ultracold  ions 
confined  in  a  ring-shaped  trapping  potential  in  a  weak  magnetic 
field.  The  mass  and  charge  of  each  ion  are  M  and  q,  respectively. 
The  diameter  of  the  ion  ring  is  d,  and  the  magnetic  field  is  B.  (c). 
Examples  of  the  pseudopotentials  (Text)  for  a  9Be+  ion  in  a 
quadrupole  ring  trap  (solid  curve)  and  a  linear  octupole  trap 
(dashed  curve)  along  the  x  or  y  axis. 


field.  The  ring-shaped  trapping  potential  can  be  created  by 
a  quadrupole  storage  ring  trap  [17,18],  a  linear  rf  multipole 
trap  [19,20],  a  multiple  trap  from  planar  ring  electrodes 
[21]  or  other  methods.  As  an  example,  we  consider  N 
identical  ions  of  mass  M  and  charge  q  in  a  ring  trap  and 
a  uniform  magnetic  field  B  [Fig.  1(b)].  The  magnetic  field 
is  parallel  to  the  axis  of  the  trap.  It  is  very  weak  so  that  it 
does  not  affect  trapping.  The  equilibrium  diameter  of  the 
ion  ring  is  d.  Figure  1(c)  shows  examples  of  the  trapping 
potentials  for  a  9Be+  ion  in  the  radial  plane  of  a  quadrupole 
ring  trap  and  a  linear  octupole  trap.  See  the  Supplemental 
Material  [22]  for  more  information. 

When  the  average  kinetic  energy  of  ions  (kBT / 2,  where 
kB  is  the  Boltzmann  constant  and  T  is  the  temperature)  is 
much  smaller  than  the  typical  Coulomb  potential  energy 
between  ions,  i.e.,  T  <C  Nq2 /{2v2eX)kBd),  the  ions  form  a 
Wigner  ring  crystal.  For  ions  in  a  ring  crystal,  we  can 


expand  the  Coulomb  potential  around  equilibrium  posi¬ 
tions  to  the  second  order.  So  the  many-body  Hamiltonian 
[22]  becomes  quadratic.  We  can  diagonalize  the  quadratic 
Hamiltonian  by  introducing  a  set  of  N  normal  coordinates 
qj  and  normal  momenta  p'- .  The  normal  coordinate  and 
momentum  of  the  collective  rotation  mode  are  qx  = 
and  p\  =  -CjLp,-,  respectively.  The  remaining 

N  —  1  normal  coordinates  correspond  to  relative  vibration 
modes.  Choosing  the  potential  energy  at  equilibrium  posi¬ 
tions  as  the  origin  of  energy,  the  Hamiltonian  of  the  system 
becomes  [23]: 


H  = 


Md2l\  dqx 


^)2+i(-£+^)} 


(1) 


where  h  =  h/(2v)  and  h  is  the  Planck  constant,  a  = 
qTrd2B/{Ah)  is  the  normalized  magnetic  flux,  r/2  = 
q2Md/ (87r/r  en),  and  oj  x  (j  >  2)  is  the  normalized  normal 
mode  frequency. 

The  lowest  normalized  relative  vibration  frequency  is 
c o2  =  2.48  when  N  =  10  and  will  increases  as  N  increases 
[23].  co2  =  V0.32AUn(0.77AO  for  large  N.  r)  =  5.6  X  104 
for  9Be+  ions  in  a  d  =  20  /xm  ring  trap.  So  it  costs  a  lot  of 
energy  to  excite  the  relative  vibration  modes.  Thus  we  have 
rij  =  0  for  all  j  >  2  at  lowest  energy  states,  where  are 
the  occupation  numbers  of  the  relative  vibration  modes. 
The  wave  function  has  to  be  symmetric  with  respect  to  the 
exchange  of  two  identical  bosonic  ions  (e.g.,  9Be+  ions), 
and  has  to  be  antisymmetric  with  respect  to  the  exchange 
of  two  identical  fermionic  ions  (e.g.,  24Mg+  ions).  For  an 
ion  ring  of  identical  bosonic  ions,  the  energy  Eni  and  the 
angular  frequency  cont  of  the  nrth  eigenstate  of  the  col¬ 
lective  rotation  mode  are  [22]: 

2Nh2 

Eni=E*(nl-a)2  =  1^(nl-a)2, 

4  h 

=  0)*(n  1  -  a)  =  (nx  -  a),  (2) 

where  E*  =  2Nh2/(Md2)  and  co*  =  4ti/(Md2)  are  the 
characteristic  energy  and  the  characteristic  frequency  of 
the  collective  rotation,  respectively.  For  identical  fermionic 
ions,  the  results  are  the  same  as  Eq.  (2)  if  N  is  an  odd 
number,  and  n x  should  be  changed  to  n ,  +  !  if  N  is  an  even 
number. 

In  classical  mechanics,  the  angular  velocity  of  the  lowest 
energy  state  is  always  w/co*  =  0,  which  means  that  the 
ions  do  not  rotate.  In  quantum  mechanics,  however, 
co/ co*  =  0  is  not  an  eigenvalue  when  the  normalized  mag¬ 
netic  flux  a  is  not  an  integer  or  half  of  an  integer.  So  the 
ions  can  rotate  persistently  even  at  the  ground  state.  Since 
the  ions  are  in  the  ground  state  already,  there  is  no  radiation 
loss  due  to  the  rotation.  The  rotation  frequency  is  indepen¬ 
dent  of  the  number  of  ions  in  the  ring.  The  energy  gap 
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between  the  ground  state  and  the  first  excited  state  is  A£  = 
Nh2 /{Md2),  when  a  =  1/4.  A£— >oo  when  V— *■  oo.  Thus 
the  persistent  rotation  of  identical  ions  is  a  macroscopic 
quantum  phenomenon  and  is  robust  for  observation  when 
N  is  large,  which  is  important  for  a  time  crystal  [3].  If  the 
relative  vibration  modes  are  not  in  their  ground  states,  the 
symmetry  requirement  of  the  c.m.  motion  of  identical 
particles  is  relaxed  and  the  energy  gap  between  different 
c.m.  motion  states  becomes  smaller.  The  result  becomes 
the  same  as  that  of  a  rigid  body  when  all  particles  are 
different  from  each  other.  Then  the  maximum  energy  gap 
between  the  ground  state  and  the  first  excited  state  is 
A£rigid  =  2ti2/(NMd2). 

Figures  2(a)  and  2(c)  show  the  lowest  energy  levels  and 
rotation  frequencies  of  an  ion  ring  consisting  identical 
bosonic  ions.  The  angular  frequency  of  the  persistent  ro¬ 
tation  of  the  ground  state  is  a  periodic  function  of  the 
magnetic  flux  [Fig.  2(c)].  Figure  2(d)  shows  the  rotation 
frequency  of  a  bosnic  ion  ring  in  a  constant  positive 
magnetic  field  fi0  as  a  function  of  its  normalized  diameter 
d/d0,  where  d{]  =  -J4h/ ( TrqB(]).  The  rotation  frequency  in 
Fig.  2(d)  is  normalized  by  co*0  =  qB0/2M.  The  ground 
state  is  n,  =0  when  d/d0<  1/ ~J2.  The  rotation  frequency 
is  independent  of  the  ring  diameter  and  the  number  of  ions 
in  the  ring  for  this  state,  and  oscillates  and  decreases  to  0 
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FIG.  2  (color  online).  The  energy  levels  and  rotation  frequen¬ 
cies  of  trapped  ions,  (a)  The  energy  levels  of  identical  bosonic 
ions  as  a  function  of  the  magnetic  flux  a.  The  quantum  number 
«!  is  labeled  on  each  curve.  The  angular  frequency  of  the 
persistent  rotation  as  a  function  a  is  shown  in  (b)  for  an  even 
number  of  fermionic  ions,  and  (c)  for  bosonic  ions,  (a)  and  (c) 
are  also  applicable  to  an  odd  number  of  fermionic  ions,  (d)  The 
angular  frequency  of  the  persistent  rotation  of  a  bosonic  ion  ring 
in  a  constant  magnetic  field  B0  >  0  as  a  function  of  its  normal¬ 
ized  diameter  d/d0. 


when  d/d0  increases  above  I  / ~J2.  If  we  confine  many  ions 
in  a  harmonic  trap  to  form  a  3D  spatial  crystal,  ions  in  the 
crystal  will  rotate  with  the  same  angular  frequency  <Uq  and 
form  a  4D  space-time  crystal  when  the  outer  diameter  of 
the  ion  crystal  is  smaller  than  d0/ -J2.  If  we  confine  ions  in 
two  concentric  ring  traps  with  diameters  larger  than 
d()/\/2,  the  rotation  frequencies  of  the  two  rings  can  be 
different  or  the  same,  depending  on  the  interaction  between 
ions  in  different  rings.  When  the  ratio  of  the  rotation 
frequencies  of  the  two  rings  is  an  irrational  number,  the 
ions  have  an  order  in  time  but  cannot  reproduce  their 
positions  simultaneously.  Thus  we  have  a  time  quasicrys¬ 
tal,  in  analog  of  a  conventional  spatial  quasicrystal  [24]. 

The  persistent  rotation  of  trapped  ions  can  be  detected 
by  measuring  the  Doppler  shift  of  moving  ions,  or  inferred 
by  probing  the  energy  levels  of  the  ion  ring.  More  impor¬ 
tantly,  we  can  observe  the  persistent  rotation  directly  by 
measuring  the  ion  positions  twice  when  N  is  large.  For 
example,  if  we  have  an  ion  ring  consisting  N  identical 
9Be+  ions  in  their  lowest  energy  state,  we  can  first  use  a 
pulse  of  two  copropagating  laser  beams  to  change  the 
hyperfine  state  of  one  (or  a  small  fraction)  of  the  ions  by 
stimulated  Raman  transition  and  use  this  ion  as  a  mark 
(qubit  memory  coherence  time  greater  than  10  s  has  been 
demonstrated  with  9Be+  ions  [25]).  Both  laser  beams  are 
parallel  to  the  axis  of  the  ion  ring  and  have  waists  of  vv0. 
We  assume  that  the  pulse  is  very  weak  so  that  on  average 
less  than  one  ion  is  marked.  This  two-photon  process 
localized  the  position  of  the  mark  ion  with  an  uncertainty 
of  about  A.r  ~  w{]/\fl.  The  amplitude  of  the  transverse 
momentum  of  each  photon  in  a  Gaussian  beam  with  waist 
of  vv(|  is  about  h/w0.  Thus  the  momentum  of  the  ion  ring  is 
changed  by  about  A  p  ~  \flfi  /  vv0 .  A  p  should  be  smaller 
than  the  absolute  value  of  the  initial  momentum  of  the  ion 
ring,  which  is  Nh/(2d)  when  a  =  1/4.  Thus  the  waists  of 
lasers  need  to  satisfy  2\fld/N  <  w0  <  \/2d  in  order  to 
localize  the  position  of  an  ion  without  significantly  alter 
the  initial  momentum  of  the  ion  ring.  This  condition  can  be 
fulfilled  when  N  is  large.  Then  we  can  use  a  global  probe 
laser  which  is  only  scattered  by  the  mark  ion  [26]  (state- 
dependent  fluorescence)  to  measure  its  angular  displace¬ 
ment  (A 0mk)  after  a  time  separation  At.  The  displacement 
of  the  mark  is  about  A 0mk  ~  «*(??!  —  a)Af  when  N  is 
large.  The  mark  ion  repeats  its  position  when  At  ~ 
27 Tl/\co*(n\  —  a)],  where  /  is  an  integer.  After  the  mea¬ 
surement,  we  can  cool  the  ions  back  to  the  ground  state  and 
repeat  the  experiment  again. 

The  process  of  observing  the  rotation  of  ions  in  the 
ground  state  also  elucidates  the  concept  of  quantum  time 
crystals  [3].  A  quantum  time  crystal  requires  the  probability 
amplitude  of  its  effective  ground  state  to  change  periodi¬ 
cally  in  time,  which  means  it  should  be  a  rotating  inhomo¬ 
geneous  state.  The  rotating  condition  is  satisfied  by  the 
none-zero  angular  momentum.  However,  for  both  the  soli- 
ton  model  [3]  and  our  ion  crystal  model,  the  probability 
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amplitude  of  the  true  ground  state  of  the  center  of  mass  is 
homogeneous,  due  to  the  cylindrical  symmetry  of  the  sys¬ 
tem.  The  required  spatial  inhomogeneity  comes  from  sym¬ 
metry  breaking,  which  can  happen  spontaneously  when  N  is 
infinite,  or  triggered  by  an  observation  (by  the  environment 
or  an  observer)  when  N  is  finite.  The  spatial  symmetry 
breaking  can  localize  the  center  of  mass  to  a  single  position 
[5]  and  makes  the  effective  ground  state  inhomogeneous 
[3].  The  resulted  rotating  inhomogeneous  state  is  not  the 
true  ground  state,  but  an  effective  ground  state  with  its 
energy  infinitely  close  to  that  of  the  true  ground  state 
when  A  — »  oo  [5],  For  both  the  soliton  model  and  our  ion 
crystal  model,  a  weak  observation  of  a  single  particle  can 
localize  the  center  of  mass  of  all  particles  and  project  the 
true  ground  state  to  a  rotating  effective  ground  state.  In 
contrast,  Cooper  pairs  or  electrons  in  a  metal  ring  have  no 
spatial  order,  so  an  observation  of  a  single  particle  can  not 
localize  the  center  of  mass  of  all  particles. 

To  study  the  effects  of  finite  temperatures  on  persistent 
rotation,  we  assume  the  trapped  ions  are  at  thermal  equi¬ 
librium  with  temperature  T.  Then  the  average  angular 
frequency  of  the  ions  is 

°°  !  I 

ft>=  £  (t^le-^/kBT)  (3) 

where  the  partition  function  is  Z  =  £„,e-£" ,/knT.  It  is 
convenient  to  define  T*  =  E* /kB  =  Nhco*/2kB  as  the 
characteristic  temperature  for  the  ring  of  ions.  T*  can  be 
considered  as  the  phase  transition  temperature  of  the  space- 
time  crystal.  It  increases  when  N  increases.  Figure  3(a) 
shows  the  average  rotation  frequency  of  an  ion  ring  con¬ 
sisting  identical  bosonic  ions  as  a  function  of  the  tempera¬ 
ture.  T*  and  of  as  a  function  of  the  diameter  of  an  ion  ring 
(or  an  electron  ring)  are  displayed  in  Fig.  3(b).  For  a  d  = 
100  fim  ion  ring  consisting  100  9Be+  ions,  T*  =  1.1  nK 
and  co*  =  2.8  rad/s.  T  is  larger  for  smaller  ion  rings.  The 
characteristic  temperature  of  an  electron  ring  is  much 
higher  because  of  smaller  mass  of  electrons. 

In  order  to  experimentally  realize  such  a  space-time 
crystal  with  trapped  ions,  we  need  to  confine  ions  tightly 
to  have  a  small  d,  and  cool  the  ions  to  a  very  low  tempera¬ 
ture.  Recently,  cylindrical  ion  traps  with  inner  radius  as 
small  as  1  /x m  have  been  fabricated  [27].  Simulations 
suggested  that  it  is  also  possible  to  confine  charged  particles 
with  a  nanoscale  rf  trap  [28],  The  challenge  is  that  ions  must 
be  cooled  to  below  1  /xK  for  a  microscale  trap  [Fig.  3(b)]. 
We  propose  to  first  add  a  pinning  potential  to  confine  ions 
with  MHz  trapping  frequencies  in  the  circumference  direc¬ 
tion.  A  combination  of  Doppler  cooling  and  resolved- 
sideband  cooling  can  be  used  to  cool  the  ions  to  the  ground 
state  of  the  MHz  trap  [15,29].  The  system  is  in  the  ground 
state  of  the  ring-shaped  trapping  potential  after  ramping 
down  the  pinning  potential  adiabatically.  For  T*  =  1.1  nK, 
the  ramping  down  time  should  be  longer  than  7  ms.  Another 
way  to  achieve  an  ultralow  temperature  of  ions  is  to  put  the 


d  (pm) 

FIG.  3  (color  online).  The  temperature  dependence  of  the 
persistent  rotation  of  trapped  ions,  (a)  The  average  angular 
frequency  of  the  persistent  rotation  of  identical  bosonic  ions  as 
a  function  of  the  temperature.  From  top  down,  the  magnetic  flux 
increases  from  —0.45  to  0.45.  (b)  The  characteristic  temperature 
(left  axis)  and  the  characteristic  frequency  (right  axis)  of  the 
persistent  rotation  of  an  ion  (or  electron)  ring  consisting  100 
identical  ions  (or  electrons)  as  a  function  of  the  diameter.  From 
top  down,  the  trapped  particles  are  electrons,  9Be+  ions,  and 
24Mg+  ions. 

ions  near  ultracold  neutral  atoms  [30],  which  has  been 
cooled  to  below  0.5  nK  by  adiabatic  decompression  [31]. 

The  ions  need  to  be  laser  cooled  to  the  ground  state  of  a 
MHz  trap  before  they  can  be  adiabatically  transferred  to 
the  ground  state  of  the  ring  trap.  The  ground  state  cooling 
of  a  MHz  trap  can  be  achieved  with  resolved-sideband 
cooling  [15,29].  Recently,  a  ground  state  population  of 
99  ±  1%  has  been  achieved  for  a  single  ion  oscillating  at 
585  kHz  in  a  room  temperature  trap  with  a  sub-Hertz 
heating  rate,  corresponding  to  a  temperature  of  only 
6  fiK  [32], 

The  electric  field  noise  on  the  trap  electrodes  should  be 
reduced  to  minimize  heating.  The  fundamental  noise  is  the 
Johnson  noise  in  the  trap  circuitry.  However,  real  experi¬ 
ments  are  dominated  by  anomalous  heating  that  is  related 
to  the  contaminations  of  the  surface  of  electrodes.  By 
cleaning  the  surface  of  electrodes  in  situ  with  an  argon- 
ion  beam,  Hite  et  al  [33]  reduced  the  heating  rate  by  2 
orders  of  magnitude  at  room  temperature.  The  heating  rate 
can  also  be  reduced  dramatically  by  reducing  the  tempera¬ 
ture  [34],  The  static  stray  electric  fields  should  also  be 
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minimized.  Recently,  the  stray  electric  field  has  been  com¬ 
pensated  to  about  0.1  V/m  for  optical  trapping  of  an  ion  in 
a  linear  trap  [35].  This  value  can  be  reduced  further  with  a 
better  design  of  the  ion  trap,  for  example,  using  only 
concentric  planar  ring  electrodes  [21].  Such  a  planar  ring 
multiple  trap  needs  only  rf  voltages  and  is  symmetric  in 
rotation. 

In  conclusion,  we  propose  a  method  to  create  and  ob¬ 
serve  a  space-time  crystal  experimentally  with  trapped 
ions.  We  also  discuss  about  how  to  create  a  time  quasi¬ 
crystal.  The  space-time  crystals  of  trapped  ions  provide  a 
new  dimension  for  studying  many-body  physics,  and  may 
have  potential  applications  in  quantum  information  for 
simulating  other  novel  states  of  matter. 
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